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§0.  Introduction  and  Summary 

Because  of  analytical  and  numerical  difficulty,  effort  is  often 
made  to  approximate  a  birth-death  process  by  a  diffusion  process.  Such 
approximation  is  usually  weakened  by  the  absence  of  either  numerical 
comparison  or  analytical  error  bound.  The  present  small  paper  was 
motivated  by  a  need  to  replace  the  Ornstein-Uhlenbeck  (0-U)  process 
by  a  certain  finite  birth-death  process,  the  Ehrenfest  process  as  des¬ 
cribed  below.  For  these  two  processes,  a  numerical  comparison  is 
possible.  On  the  one  hand,  extensive  numerical  results  for  the  0-U 
process  are  available  [8,9].  Simultaneously,  algorithms  and  software 
available  for  birth-death  processes  [5,7]  extend,  with  some  further 
analytical  and  numerical  effort,  to  the  approximating  Ehrenfest  process. 

The  Ehrenfest  process,  a  finite  Markov  chain  of  some  practical 

interest,  arises  from  the  sum  of  K  independent  identical  chains  Jj(t), 

each  on  state  space  {0,1}  and  governed  by  transition  rates  v„,  =  v 

K 

1/2,  The  Markov  chain  N„(t)  =  J  J.(t)  has  transition  rates  v 


1 


lin  N^(t)  =  J  has  transitic 

j  =  ll  ^ 


n,n+l 


y(K-n),  0  <  n  <  K-1,  and  =  y  1  s  n  <  K.  Consequently,  the 

K 

local  growth  rate  of  the  variance  v  ^  - 

°  n,n+i.  n,n-l  ^ 

^v(Nj.g(t) ,  Nj.^(t+T)]  = 


KS' 


associated  stationary  chain  Nj^g(t),  ore 

K  ”  T 

j  e  and  asymptotic  normality. 

These  simple  properties  suggest  that  the  sequence  of  processes 
XyCt)  =  /I  N2v(t)  -  converges  in  law  to  the  Ornstein-Uhlenbeck 
process  as  V  -♦  *,  since  such  0-U  processes  are  characterized  by  their 
Markov  property,  normal  distribution  and  exponential  covariance  function. 
The  object  of  this  research  note  is  to  establish  this  convergence  in  law. 
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to  develop  systematically  the  properties  of  N2y(t),  and  to  quantify 
its  behavior. 

Of  particular  interest  is  the  evaluation  of  the  quality  of  the 

chain  X^Ct)  as  an  approximation  to  the  0-U  process.  Even  though  the 

0-U  process  X^(t)  is  quite  tractable  [8,9],  simple  related  processes 

such  as  Yg(t)  =  6  /  e  ^  ^X^(t')dt'  are  intractable.  The  latter 
-00 

process  Y0(t)  describes  an  averaged  0-U  process  with  averaging  time  6 
Evaluation  analytically  of  the  passage  time  distributions  for  Y^Ct)  and 
related  distribution  of  the  maximum  of  Y.(t)  over  an  interval  has  defeated 
effort  to  date. 

Recently  the  authors  have  developed  an  algorithmic  framework  [10] 
for  quantifying  the  behavior  of  finite  bivariate  Markov  chains  [Nj(t), 
N^Ct)]  where  the  marginal  process  N2(t)  is  lattice-continuous.  The 
bivariate  Markov  process  [X  (t),  Y  (t)]  may  be  approximated  by  a  finite 
bivariate  Markov  chain  [X  (t),  Y  (t)]  whose  sequence  converges  in  law 
to  [X^(tJ,  Yg(t)].  It  is  expected  that  the  use  of  Xy(t)  in  conjunction 
with  that  framework  will  permit  evaluation  of  the  distributions  needed 
for  Y(t) . 

The  structure  of  this  paper  is  as  follows.  Following  Karlin  and 
McGregor  [2],  notation  is  established  for  the  spectral  representation  of 
birth-death  processes  in  §1.  The  specific  form  appropriate  to  the 
Ehrenfest  chain  N2yCt)  is  developed  in  §2.  Related  first  passage  time 
density  structure  is  described  in  §3.  In  §4,  the  first  passage  time 
structure  for  the  Ehrenfest  chain  N2y(t)  is  studied  in  detail.  In  par¬ 
ticular,  for  the  first  passage  time  T^y  of  N2y(t)  from  0  to  V,  we  show 


4 


that  21^^  -  log  V  converges  in  law  to  a  distribution  conjugate  to  the 
extreme-value  distribution.  The  convergence  in  law  of  Xy(t)  to  the  0-U 
process  and  related  results  are  established  in  §5.  Finally,  in  §6,  the 
numerical  results  are  exhibited.  In  particular,  the  ergodic  probabili¬ 
ties  and  the  first  passage  time  distributions  of  Xy(t)  for  V  =  10,  25 
and  40  and  those  of  the  0-U  process  are  compared. 


§  1 .  The  spectral  representation  of  birth-death  processes 

A  birth-death  process  N(t)  is  a  Markov  chain  in  continuous  time 
on  n  =  {0,1,2,...}  governed  by  a  transition  rate  matrix  y  where 


X  >0  n=m+l,  m=0,l,2,.. 
m  ’  >  >  » 


M^>0  n=m-l,  m=l,2,, 


n-m  >  1 


Let  Q  =  V  -  Vp  where  is  a  diagonal  matrix  with  m-th  diagonal  element 

equal  to  J  v  .  The  transition  probability  matrix’  P(t)  is  then  given 
ntN  ™ 

by  [5] 


(1.1)  P(t)  =  e‘ 


where  I.  The  matrix  Q  is  called  the  infinitesimal  generator  of  the 

chain  N(t) . 

In  a  scries  of  papers  [2],  [3],  [4],  Karlin  and  McGregor  analyze  the 
spectral  representation  of  P(t)  and  use  the  results  to  study  various 
probabili.stic  quantities.  In  this  section  we  briefly  describe  their 
results  and  establish  notation. 

The  infinitesimal  generator  Q  has  a  vector  eigenfunction  y(x)  = 
with  associated  eigenvalue  -x,  i.e.. 


(1.2)  Qy(x)  =  -xy(x)  , 
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or  componentwise  one  has  with  YqCx)  =  1, 


(1.3) 


We  note  from  (1.3)  that  y  (x)  is  a  polynomial  of  degree  n  with  leading 

- 1 

coefficient  (11  X.)  .  Let  f(x,t)  be  a  vector  function  defined  by 

j=0  J 


(1.4)  f(x.t)  =  P(t)y(x)  . 

From  the  Kolmogorov  forward  equations,  one  then  has 


f(x.t)  =  ^  P(t)-y(x)  »  P(t)Qy(x)  =  -xP(t)y(x)  , 

i.e. , 

(1.5)  1^  f(x,t)  =  -xf(x,t)  . 

Since  f(x,  0+)  =  y(x) ,  the  partial  differential  equation  (1.5)  is  solved 
by 

(1.6)  f(x,t)  =  e"’‘^y(x)  . 

Equations  (1.4)  and  (1.6)  imply  that  f(x.t)  =  P(t)y(x)  =  e'*^y(x).  i.e., 

K  Pmn^'^yn^^)  =  m  e  N  . 


(1.7) 


It  is  known  [3]  that  there  exists  a  measure  on  [0,®)  such  that 

tyjj(x)  becomes  a  set  of  orthogonal  polynomials  with  respect  to  tKx) . 
One  has 


Cl. 8)  /  y  Cx)y  (x)d4'Cx)  =  6  /n  ,  m,n  e  N 

where  are  the  potential  coefficients  [15]  with  =  1  and  it  = 
n-1  "n  " 

n  A./  n  p..  From  (1.7),  {p  (t) }  ^  may  be  recognized  as  the  general 

J =0  ^  J“1  ^  ncN 

ized  Fourier  coefficients  of  f^(x,t)  associated  with  the  orthogonal  poly 

nomials  {y  (x) }  for  each  m  e  W.  One  then  has,  from  (1.7), 
n  nc/v  V  j 

00 

Pmn^^^  '  0  ’  m,n  e  W  . 
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§ 2 .  The  spectral  representation  of  the  Ehrenfest  process 

Consider  2V  independent  and  identical  Markov  chains  J^(t)  in  con¬ 
tinuous  time  on  {0,1},  governed  by  the  transition  ra<-e  matrix 

0^'  transition  probability  matrix  of  (t)  . 

One  finds  that 


(2.1)  9(t)  = 


.o' 

o 

o 

qoi(t) 

f(t) 

g(t) 

Lqio^^^ 

qil(t)_ 

_g(t) 

f(t)_ 

where  f(t)  =  -^(1  +  e  ^)  and  g(t)  =  ^{1  -  e  ^) .  For  future  reference  we 
define  the  generating  functions 

(2.2)  aQ(t,u)  =  qQo^^) 

and 


(2.3)  ai(t,u)  =  qii(t)u  =  q(t)  +  f{t)u 


11' 


2V 


Let  N2y(t)  =  J  >l^(t).  Then  N^y(t)  is  a  birth-death  process 


j  =  l 


1 


2V' 


on 


N  =  {0,1, . . . ,2V} ,  governed  by  the  transition  rates 

(2.4)  =  i(2V-n)  ;  t  «  N  . 


We  note  that 


(2.5) 


'^n  =  ^n  "  %  =  ^  ^  • 


This  birth-death  process  may  be  called  the  Ehrenfest  process  [1].  Let 


P(t)  be  the  transition  probability  matrix  of  N2y(t}  and  define  the  gen¬ 
erating  function 


2V 

(2.6)  ^  ,  m  e  W  , 

k=0 

From  the  independence  of  the  (t) ,  one  has 

(2.7)  6jj,(t.u)  =  aQ(t.u)^''''"aj(t.u)'"  =  {f(t)  +  g(t)u}  "^'''"{g(t)  +  f(t)u}"' 

It  can  be  readily  seen  from  (2.7)  that  p  (t)  is  a  linear  combination  of 

mn 

“  kt 

{e  ),  0  s  k  <  2V.  Karlin  and  McGregor  [  ]  identify  the  orthogonal  poly¬ 
nomials  and  associated  discrete  measure  described  in  §1  for  the  Ehrenfest 
process  and  give  the  following  spectral  representation. 

1  2V 

(2.8)  y  (X)  =  I  (  „'^)(^)(-l)^  .new. 

(;; )  j=0  "'J  ^ 

(2.9)  dij/(x)  =  ,  x  =  0.1,. ..,2V  , 

and 

2V 

2V  -2V  r  2V  -it 

(2.10)  P  (t)  =  (  )2  y  C^)y  (i)y  (j)e  . 

'  ^  mn'  '  'n  i  -^m 

j=0 

The  polynomials  (yj^(x))  are  called  the  Krawtchouk  polynomials. 

The  ergodic  distribution  of  Jj(t)  is  (-j,  ^) .  Therefore  the  ergodic 
T 

distribution  e  =  [eQ,...,e2y]  of  the  Ehrenfest  process  is  given  by 


C2.ll) 

from  the 

(2.12) 

as  t  ->■  “ 


.2V,.-2V  ,, 

en  =  Cn  )2  .  n  c  W 


independence  of  J^(t).  We  notice  that 


2V 

^2V,--2V  r  .2V,  ...  -it 

p  (t)  =  (  )2  )  (•  )y  (i)v  (i)e  -*■  e 

j=0  •’ 


as  expected. 
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§3.  The  first  passage  time  and  the  historical  maximum  of  birth-death 
processes 

Let  be  a  birth-death  process  on  n  =  {0,1,..., K}  (K  <  ®) . 

Define  the  first  passage  time  T  of  N„(t)  from  state  m  to  state  n 
by 

(5.1)  T  =  inf{t':  N,.(t')  =  niN(O)  =  m)  . 

mn  K  ivy 

Let  s  (t)  be  the  p.d.f.  of  T  with  Laplace  transform  o  (s) .  For 
mn  mn  mn 

notational  convenience  we  denote  T  ,  by  and  correspondingly  s^(t) 

m,m+l  m  ^  '’  m 

and  a*(s) .  From  a  probabilistic  argument  [5],  one  has  the  following 
consistency  relations 

(3.2)  o^(s)  =  X^[s  +  -  %oj^_j(s)]‘^  ,  n  =  1,2,.. .,K 

Ay 

wit})  a.(s)  =  — ,  Let  {g  (s)}  be  the  polynomials  of  order  n  defined 

u  s  ^  Aq  n 

by 


(5.3) 


(s)  =  — r~^ 
On'  g^(s) 


,  n  =  1,2, ... ,K  ;  gy(s)  =  1 


Since  o„  (s)  =  o,,  ,  (s)o  ,  (s) ,  we  have  from  (3.2) 

On  On-l  n-1  ’ 


^  ^  0,1,2,.. .,K  , 

m 

where  g  ^^(s)  =  0.  As  shown  in  (  ],  the  polynomials  gj^(x)  are  related 

to  the  orthogonal  polynomials  y^^Cx)  given  in  §1  by 

(3.5)  y^ix)  =  gj,(-x)  . 
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We  notice,  from  (3.3),  that 


(3.6)  0^(5)  =  g^(s)/g^^j(s)  . 

From  (5.5),  {gj^(s)}  are  orthogonal  polynomials  anJ  hence  the  zeros  of 
g^(s)  are  distinct,  and  the  zeros  of  any  two  successive  polynomials 
interleave.  Furthermore,  the  zeros  are  negative.  Therefore  is 

of  the  form 


(5.7)  o*(s)  =  J  p  . 

n  ^ni  s  +  o  . 

j=0  ■>  nj 


i  1 

where  p^j  >0,  J  p^^  =  1  and  are  the  zeros  of  gp+j(s).  Hence 


j=0 


Sjj(t)  is  completely  monotone.  A  more  detailed  discussion  can  be  found 
in  Keilson  [  ].  The  downward  passage  time  n  1  ”  ^n  treated 

similarly  when  K  <  ®. 

Let  M|^(nQ,  6)  be  the  historical  maximum  of  N|^(t)  in  the  interval 
[0,6]  given  that  N|^(0)  =  n^,  i.e.. 


(3.8)  M  (n  ,  e)  =  max  (N  (t)iN..(0)  =  n.) 

^  ^  o<t<e  ^  0 


We  observe  the  following  dual  relation  between  the  first  passage  time 

'"O  ‘  “"■*  ”k'"o-  *>• 


15 


Hence  one  has  for  integers  n 


(3.9)  F  g(n)  = 


For  the  corresponding  stationary  process  c.d.f.  of  the 

historical  maximum  is  given  by 

(5.10,  F,(„)  =  I  .5  (6) 

m<n 

T 

where  e  =  (e  ]  is  the  stationary  distribution  of  N'  _(t). 

“  n  Kd 
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S4.  The  first  passage  time  structure  of  the  Ehrenfest  process 

As  we  saw  in  §2,  the  Ehrenfest  process  N^yCt)  is  a  birth-death 
process  on  {0,1,. ..,2V}  governed  by  transition  rates  =  (2V-n/2, 

0  <  n  s  2V-1,  and  =  n/2,  1  s  n  5  2V.  The  recursion  (3.4}  then 
becomes 


(4.1) 


2V 


—  [(sn')g^(s)  -  ^  g^.iCs)],  0  <  n  <  2V 


with  g  j(s}  =  0  and  gQ(s}  =  1.  From  (2.8)  and  (3.5),  gj^(s)  are  given 
explicitly  by 

2V 

(4.2)  ^  i  .  0  s  n  <  2V  . 

"  (  )  n=0  ^ 

In  order  to  evaluate  the  first  passage  times  s  (t)  (m<n)  with 

Laplace  transform  o^(s)  =  tJ^(s) . °n-l^^^  "  Pieros 

of  gjj(s)  are  needed.  The  amount  of  effort  required  for  the  zero  search 
is  considerably  reduced  by  the  following  observations. 

Theorem  4.1 

Let  h^(s)  =  gj^(s-V).  Then  hj^(s)  =  (-l)"h^(-s),  n  >  0,  i.e., 

{odd  when  n  is  odd, 
even  when  n  is  even. 


Proof 

Eq.  (4.1)  can  be  written  in  terms  of  h^^Cs)  by 


15 


i 


and  the  result  follows  by  induction  on  n.  C 
Corollary  4,2 

(a)  g„(-x)  =  0  =>  g  (x  -  2V)  =  0 

ll  It 

(b)  n:  odd  =>  g^(-V)  =  0 

Proof 

This  follows  immediately  from  Theorem  4.1.  D 
Theorem  4.1  implies  that  the  zeros  of  h^(s)  are  symmetric  about 
zero  and  correspondingly,  as  in  Corollary  4.2,  the  zeros  of  are 

symmetric  about  -V.  Hence  we  need  to  find  only  [(n-l)/2]  zeros,  where 
[x]  is  the  minimum  integer  which  is  greater  than  or  equal  to  x.  Further¬ 
more,  there  are  only  1  +  [(n-l)/2]  terms  in  each  h^(s)  while  gjj(s)  has 
(n+1)  terms.  We  therefore  reduce  computation  time  approximately  by  a 
factor  of  4.  We  note  that  similar  results  are  available  for  general 

birth-death  processes  whenever  v  =X  =vis  constant. 

n  n  n 

We  next  show  that  gy(s)  has  negative  odd  integers  as  its  roots. 

A  preliminary  lemma  is  needed. 

Lemma  4 . 3 

gn(''"J  =  Sjn(-n)  ,  0  <  m,n  <  2V,  m,n  integers  . 

Proof 

This  is  a  known  result  [  ].  A  simple  proof  is  as  follows. 
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2V 


(m-j) !  (2V-n-ni+j)  !  j!(n-j)! 

2V-»,  . 


(2V) ! 


2V 


and  the  result  follows  from  (4.2).  □ 

Theorem  4 . 4 

V 

g  (s)  =  n  (s+2j-l)/(2j-l) 

j  =  l 


Proof 

Corollary  4,2  (b)  states  that  g^^C-V)  =  0  for  n  odd.  The  theorem 
is  then  immediate  from  Lemma  4.3.  0 

Eq,  (3.3)  along  with  Theorem  4.4  enables  one  to  write  the  Laplace 
transform  of  the  first  passage  time  T^^  from  0  to  V. 


(4.3) 


V 

n 

j=i 


s+2j-l 


s+2j-l 


where  c . 
J 


n 

k-1 

k#j 


,2V 


1 

‘  2j-l 


The  limiting  distribution  of  T^^  can  now  be  evaluated.  We  refer  the 
reader  to  Keilson  [6]  for  the  definition  of  the  conjugate  distribution. 
Theorem  4.5 

2TQy  -  log  V  ^  G  as  V  +®  where  G  corresponds  to  the  distribution 
conjugate  to  the  extreme-value  distribution  with  p.d.f. 
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p.  r  r ) 


/„ 


o 


-e 


-T 


_  «•  -y  <00 


Proof 


The  p.d.f.  of  2Tq^.  -  lo^  V  has  the  Laplace  transform 


V^o  (2sJ  =  V®  11  — 

OV  2s+2j-l 


By  simple  algebra  this  can  be  rewritten  as 


(4.4) 


V^o  (2s)  =  - —  n  ^ } 

Ov'-  -2V  '• 

2  3  =  1  •' 


The  factor  inside  braces  converges  to  r(s+h)  as  V  +»,  while  the  rest 
converges  to  IZ/tT  =  l/r(i5),  i.e.. 


It  is  known  that  r(l+s)  is  the  Laplace  transform  of  the  extreme  distribu- 

-t 

tion  with  p.d.f.  e  e  ^  -00  <  <  ao^  and  thus  the  theorem  follows.  □ 

We  saw  in  (3.7)  that  the  upward  passage  time  T*  is  a  finite  mixture 

of  exponential  variates.  The  fact  that  the  exit  rate  of  each  state  is 

constant,  i.e.,  v  =  V,  n  e  N,  enables  one  to  show  that  T*  can  also  be 
n  n 

expressed  as  an  infinite  mixture  of  Gamma  variates  of  odd  order. 

Theorem  4.6 

T*  and  T^  for  the  Lhrenfest  process  are  infinite  mixtures  of  Gamma 
variates  of  odd  order,  r(V,  Zj+l),  j  =  0,1,2,...  . 
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Proof 


The  recurrence  formula  for  o^{s)  in  (3.2)  can  be  rewritten  as 


(4.5) 


^  -  ^n 


c(s)  =  Oq(s) 


V 

s+V 


where  r^  =  1  -  ^  is  the  probability  of  going  up  given  exit  from  n 
and  r”  =  ^  is  that  of  going  down.  For  Re(s)  >  0,  Eq.  (4.5)  has  a 
series  expansion 


"  j=0  " 

and  the  result  follows  by  induction  for  o^(s).  For  ojj(s)  it  suffices 

to  note  that  o*(s)  =  ol,,  (s) .  □ 

n  jv-n 


following  theorem  we  prove  that  when  N^ytO)  is  chosen  appropriately 
Xy(t)  converges  in  distribution  to  for  every  fixed  t,  t  2  0, 
as  V  -*•  “>.  For  notational  convenience  define 


(5.5) 


ny(x)  =  [ 


x] 


where  [x]  is  the  minimum  integer  which  is  greater  than  or  equal  to  x. 
Theorem  5 . 1 

Let  X^yCt)  be  the  0-U  process  with  X^^CO)  =  x^,  -“  <  x^  <  Let 

fl 

X^(t)  be  as  in  (5.3)  with  Xy(0)  =  y  ^  where  V  is  chosen  large 

enough  so  that  -/2V  <  X^(0)  <  .  Then  Xy(t)  ^  XQu(t),  for  all  t, 

t  2  0,  as  V  +». 


Proof 

-wX  (t)  /y 

Let  ()iy(x^,w,t)  =  E[e  1X^(0)  =  One  sees  from 

(2.6)  and  (5.3)  that 


(5.6) 


4ly(XQ,W,t) 


w/2V 

e 


where  N2y(0)  =  V  +  ny(xQ).  We  wish  to  show  that  (fiy(Xp,w,t)  -»•  Y(Xp,w,t) 
as  V  Eq.  (5.6)  can  be  rewritten,  by  (2.7),  as 


w/zv 

<l'y(xQ,w,t)  =  e  [(f(t)  +  g(t)e 


■^/v 


)(g(t)  +  f(t)e 


V 


g(t)  *  f(t)e 


/2 

-w/y 


nv(xQ) 


f(t)  +  g(t)e 
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Since  f(t)  +  g(t3  =  1,  the  first  factor  in  the  above  equation  becomes 


e'"'^6y(t,  e  =  (1  +  2f(t)g(t){cosh(w/|)  -  1}]^ 


/2  1 
u)  ■  1 1  ^  <  •5-  so  that 


logfe'^'^BvCt,  e  =  V  log[l  ♦  2f (t)gCt) {cosh (w /|)  -  D] 


V  I  i(2f{t)g(t))^{co5h(w/|)  -  1}^ 


k=l 


1  ,,  -2t,  2 

y  (1  -  e  )w  +  0(V  ) 


The  second  factor  can  also  be  rewritten  as 


-/I 

1  .  f(t)  -  g(t)  (1  -  e  ") 


-w 


L 


f(t)  +  g(t)e 


r  e'^^w  +  0(V'^) 

1 - - 

1  +  0(V 


\^V 


1  -  - 


x^c*^w  +  XqO(V' 


x^yi  il  ^  0(V'^)} 


0»'  2 


'Or 


/I 
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Therefore  the  second  factor  converges  to  exp{-  x^e  )  as  V  >■  o".  As 
a  result,  one  has 

-t  J  -2t 

({'^,U,^,w,t)  ->  exp{-  XjjC  w  +  -  (1  -  e  )w  } 
as  V  ”,  proving  the  theorem.  D 

The  follov>ing  corollary  is  immediate  from  Theorem  5.1. 

Corollary  5.2 

For  gi\-cn  x^,,  x  real,  let  m  =  V  +  n,,(x,j  and  n  =  V  n  .(x)  .  Then 
^  C'  V  U  V 


p__ ( t ) 
mn 


g(Xp,x,t )  as  V  O' 


for  even-  t,  t  >0. 

Corollary  5.2  ma)  be  seen  alternatively  in  the  following  manner.  K’c  first 
note  that  gfxQ.x.t)  in  (5.1)  has  a  scries  representat ion  |  ] 


(5.7) 


g(x,,x,f)  = 


illl. 

' 


w 

).  He  (x„)He  (x)  iy  c 

j=0  J  i'- 


-it 


where  {H,,  .(x)}”  is  the  set  of  the  standardized  Hermite  polynomials.  From 
^1  '  C' 

(2.10)  and  Lemma  4.5,  one  obtains  that 


(5.8) 


2V 

/V  /V  ^  n^ 

2 


,2V 


y  (^^)y.(ra)y.(n)e’^^ 
j=0  •'  ■'  ^ 


\  1 

It  is  known  (Szego  [11],  p.  -'^6)  that 


lim  /y2V  >-.(n,  =  — —  He.(x  )  , 

v-H»  /  i  '  t'TT 


and  similarl)- 


lim  /(  J  yJn) 

V-w  •'  •* 


He^ (X)  . 


The  first  factor  in  (5.8)  converges  to  e  as  V  ->■  “ 

/v 

from  Stirling's  formula  and  gC^Q.x.t)  as  V 

It  is  natural  to  expect  that  the  first  passage  time  distribution 
of  X^.(t)  also  converges  in  distribution  to  that  of  XQjj(t).  We  next  prove 
this  formally. 

Theorem  5.5 

Let  m  and  n  be  as  in  Lemma  5.2.  Let  T  ,  .  ,  .  be  the  first  passage 

r(m)r(n)  ^  ^ 

time  of  X^^(t)  from  r(m}  to  r(n)  and  let  T^  ^  be  that  of  X^^jCt)  from  Xq 
to  X.  Then 


T  .  ^  ^  T  as  V  ^  ®  . 

r(m)r(n)  x^x 


Proof 


The  Laplace  transform  of  the  p.d.f.  of  T  is  given  [  ]  by 

XpX 


-sT  r 

x„x  y(x-,x,s) 

1  ■  ffxTiT^r 


It  is  clear  that  T  ,  .  ,  >  =  T  where  T  is  the  first  passage  time  of 
r(m)r(n)  mn  mn  r  & 

N-,,(t)  from  m  to  n.  Let  s  (t)  be  the  p.d.f.  of  T  with  Laplace  trans- 
2V  mn  mn 

form  0  (sj.  From  consistency  relations,  one  has 
mn  ^  ' 


(5.10)  s  (t)*p  (t)  =  p  (t) 

'■  ^  mn  ^  ^  *  nn  ^  ^  *  mn  '• 


§f).  Nuiiicrical  Ucsu]t_s 

In  this  section  we  compare  some  results  for  the  Ehrenfest  process 
with  those  obtained  for  the  0-U  process.  The  results  for  the  0-U  pro¬ 
cess  are  described  in  [8J.  Those  for  the  Ehrenfest  process  were  derived 
in  a  two-step  procedure.  First,  the  analysis  of  the  earlier  sections 
was  used  to  obtain  numerical  answers  at  the  mass  points  of  the  Ehrenfest 
distribution.  Then,  Lagrange  interpolation  was  used  to  extrapolate  to 
the  actual  points  desired. 

Table  6.1  displays  the  ergodic  cumulative  distribution  function 
Fy(x)  =  s  x]  for  V  =  10,  25,  40,  and  80.  The  last  column  gives 

X  2  . 

that  of  X^ij(t),  i.e.,  F(x)  =  P[Xqjj(“)  s  x]  =  /  {e”^^  •  - }dy.  The 

interpolation  procedure  used  here  is  as  follows.  For  a  given  real  x, 
we  find  an  interval  [r(m),  r(m+l))  so  that  r(m)  s  x  <  r(m+l).  Then 
Lagrange  interpolation  at  the  points  {r(m-l),  r(mj,  r(m+l),  r(m+2)} 
is  used.  As  expected  from  the  central  limit  theorem,  the  error  here 
is  quite  small,  less  than  10  ^  for  V  =  80. 

In  tables  6.2(a)  and  (b) ,  the  mean  and  standard  deviation  of  the 


upward  first  passage  time  T  and  the  downward  first  passage  time  T..„ 
are  shown.  For  each  x,  the  corresponding  row  shows  the  mean  of  T^^^ 

(or  T^^)  for  X^,(t)  with  V  =  10,  25,  40,  80  and  for  X^^Ct) .  The  second 
row  shows  the  associated  standard  deviation.  For  X^,(t) ,  the  recurrence 
formulas  obtained  from  (3.2)  (cf.  Chapter  5  of  [5])  are  used  and  then 
an  interpolation  procedure  as  above.  For  the  mean  first  passage  time 


E[Tox]>  V'  =  40  gives  3  digit  accuracy  for  .1  <  x  £  2.2,  while  V  =  80 
provides  the  same  accuracy  for  .1  <  x  <  2.6.  The  relative  error  increases 
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as  X  grows  larger  than  3,  probably  due  to  numerical  imprecision. 

For  the  mean  downward.s  first  passage  time  both  V  =  40  and 

V  =  80  provide  3  digit  accuracy  for  .1  s  x  <  7.1. 

The  remaining  tables  6.3(a)-(c)  compare  the  cnnulative  distri¬ 
bution  functions  of  the  first  passage  times  for  V  =  10,  2.S.  40  and 

for  XQ^(t).  Except  for  small  values  of  t  and  small  cumulative 

probabilities,  agreement  is  quite  good,  even  at  V  =  25.  For  example, 
the  relative  error  for  jCt),  t  >  .75  is  less  than  O.l®  (for  V  =  40). 

For  V  =  25  the  same  error  is  less  than  0.15-o.  For 

bers  are  1.7'<.  for  V  =  40  and  2.1%  for  V  =  25. 

The  greatest  difficulty  with  the  numerical  procedures,  described 
in  the  previous  sections,  occurs  when  attempting  to  find  the  zeros 
of  the  orthogonal  polynomials.  Of  course,  the  major  problem  is  that, 
even  with  the  simplification  allowed  by  Corollary  4.2,  the  p""  j '“.ami.- ’ 
ggQ(x)  is  of  degree  40.  The  recursive  definition  of  (3.4  may  be  u.sed, 
at  somewhat  greater  expense,  as  a  more  stable  method  of  evaluation. 
Instability  using  a  straightforward  polynomial  evaluation  is  detected 
by  a  failure  to  find  all  zeroes,  thus  indicating  that  the  zero  bound¬ 
aries,  i.e.,  the  previous  zeroes,  are  incorrect. 

Once  the  zeroes  arc  all  found  (a  procedure  which  need  be  done  only 
once),  evaluation  of  passage  time  and  exit  time  distributions  is  straight¬ 
forward  and  relatively  inexpensive. 

Note:  The  missing  data  in  table  6.2  is  due  to  two  factors.  First, 
the  state  .spacc.s  do  not  cover  the  range  comjjletely  for  V  -  10,  25,  and 
40.  Second,  some  variance  calculations  exceeded  the  10^^  maximum  real 
number  on  the  computer. 
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X  10  25  40  80 


0.00000 

0.20000 

0.40000 

0.60000 

0.80000 

1.00000 

1.20000 

1.40000 

1.60000 

1.80000 

2.00000 

2.20000 

2.50000 

2.80000 
3.00000 
3.20000 

3.50000 

3.80000 
4.00000 
4.30000 
4.60000 
5.00000 
6.00000 


0.50000 

0.57772 

0.65149 

0.72102 

0.78239 

0.83564 

0.87971 

0.91476 

0.94179 

0.96174 

0.97565 

0.98534 

0.99352 

0.99745 

0.99873 

0.99941 

0.89983 

0.99995 

0.99999 

1.00000 

1.00000 

1.00000 

1.00000 


0.50000 

0.57842 

0.65389 

0.72363 

0.78575 

0.83901 

0.88281 

0.91748 

0.94387 

0.96314 

0.97669 

0.98579 

0.99373 

0.99747 

0.99869 

0.99935 

0.99979 

0.99994 

0.99998 

0.99999 

1.00000 

1.00000 

1.00000 


0.50000 

0.57870 

0.65440 

0.72441 

0.78666 

0.83987 

0.88360 

0.91815 

0.94437 

0.96351 

0.97691 

0.98591 

0.99375 

0.99746 

0.99867 

0.99934 

0.99978 

0.99994 

0.99597 

0.99999 

1.00000 

1.00000 

1.00000 


0.50000 

0.57899 

0.65491 

0.72506 

0.78738 

0.84060 

0.88426 

0.91869 

0.94479 

0.96379 

0.97708 

0.98601 

0.99377 

0.99745 

0.99866 

0.99933 

0.99978 

0.99993 

0.99997 

0.99999 

1.00000 

1.00000 

1.00000 


0.50000 

0.57847 

0.65415 

0.72445 

0.78719 

0.84088 

0.86489 

0.91944 

0.94544 

0.96424 

0.97733 

0.98610 

0.99377 

0.99743 

0.99865 

0.99932 

0.99977 

0.99993 

0.99997 

1.00000 

1.00000 

1.00000 

1.00000 


Table  6.1.  The  ergodic  cumulative  distribution 
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X 

10 

25 

40 

80 

3.0000 

-3 

3.9390 

-3 

6.4884 

-2 

9.987E 

-2 

1 ,2217 

^  J 

1.3055 

~2 

3.2359 

0 

0.0000 

0 

£.9620 

-2 

1 .6019 

-1 

1.9065 

-I 

2.0000 

-1 

4.2376 

-1 

2.4612 

-3 

2.6674 

-1 

2.7128 

•1 

2.7248 

-1 

2.2535 

0 

3.1277 

-1 

4.0158 

-1 

4.1673 

-1 

4.2056 

-1 

3.0000 

-1 

4.9407 

-1 

4.2235 

-1 

4.2326 

“1 

4.2606 

-1 

4.2739 

-1 

1.5926 

0 

6.8453 

-1 

6.8313 

-1 

6.9730 

-1 

6.9966 

-1 

n.oooo 

-3 

6.0283 

-1 

5.8724 

-1 

5.9280 

-1 

5.9463 

-1 

5.9720 

-1 

1.2532 

0 

9.9327 

-1 

1.0282 

0 

1.0331 

0 

1.0408 

0 

5.0000 

-3 

7.8485 

-1 

7.7414 

-1 

7.7734 

-1 

7.8128 

-1 

7.8416 

-1 

3.9231 

0 

1.4244 

0 

1.4313 

0 

1.4515 

0 

1,4608 

0 

6.0000 

-1 

1.0020 

0 

9.7759 

-1 

9.8296 

“1 

9.8754 

-1 

9.9099 

-1 

2.8483 

0 

1.9175 

0 

1.9494 

0 

1.9715 

0 

1.9819 

0 

7.0000 

-1 

3.2225 

0 

1.2012 

0 

1.2112 

0 

1 .2160 

0 

1.2205 

0 

3.4017 

c 

2.5125 

0 

2.5879 

0 

2.6109 

0 

2.6336 

0 

8.0000 

-3 

1.4582 

0 

3.4596 

0 

3.4648 

0 

1  .4729 

0 

1.4780 

0 

3.8188 

0 

3.3718 

0 

3.3802 

0 

3.4367 

0 

3.4552 

0 

9.0000 

-3 

1.7235 

0 

3.7392 

0 

1.7548 

0 

1.7599 

0 

1.7668 

0 

4.4481 

0 

4.2845 

0 

4.4523 

0  • 

4.4654 

0 

4.5004 

0 

1.0000 

0 

2.0608 

0 

2.0583 

0 

2.0730 

0 

2.0844 

0 

2.0934 

0 

6. 9835 

0 

5.5310 

0 

5.6787 

0 

5.7801 

0 

5.8420 

0 

1.2000 

0 

2.8383 

0 

2.8374 

0 

2.8,600 

0 

2.8777 

0 

2.8902 

0 

3.1026 

1 

9.1677 

0 

9.5154 

0 

9.7474 

0 

9.8567 

0 

1.3000 

0 

3.3003 

0 

3.3202 

0 

3.3528 

0 

3.3656 

0 

3.3804 

0 

1,2905 

1 

1.2050 

1 

1.2612 

1 

1.2708 

1 

1.2867 

1 

1.6000 

0 

5,3949 

0 

5.3005 

0 

5.3449 

0 

5.3772 

0 

5.3988 

0 

4.7614 

1 

2.7571 

1 

2.8596 

1 

2.9436 

1 

2.9800 

1 

2.0000 

0 

1.1064 

1 

1.0283 

1 

1.0384 

1 

1.0390 

1 

1.0428 

1 

4.3982 

2 

9.5485 

1 

1.0398 

2 

1.0348 

2 

1.0527 

2 

2.2000 

0 

1.5098 

1 

1.4755 

1 

1.4873 

1 

1.4923 

1 

1.4935 

1 

3.0310 

2 

1.9624 

2 

2.0883 

2 

2.1352 

2 

2.1448 

2 

2.6000 

0 

3.8511 

1 

3.2698 

1 

3.3363 

1 

3.3462 

1 

3.3472 

1 

7.3992 

3 

6.0950 

2 

9.8385 

2 

1.0477 

3 

1.0848 

3 

3.0000 

0 

1.6165 

2 

8.6566 

1 

8.7087 

1 

8.8338 

1 

8.6931 

1 

6.9182 

5 

3.5775 

3 

5.5132 

3 

7.6500 

3 

7.4243 

3 

3.5000 

0 

1.5634 

o 

3.5555 

2 

3.8293 

2 

3.7635 

2 

3.6436 

2 

1.0113 

5 

1.2821 

5 

1.3202 

5 

4,0000 

0 

2.3792 

3 

2.1492 

3 

2.0184 

3 

4.8987 

6 

3.5357 

6 

4.0691 

6 

5,0000 

0 

1.8053 

5 

1.7567 

5 

1.4074 

5 

1.8622 

10 

1.9808 

10 

6.0000 

0 

5.1941 

7 

2.8268 

7 

2.4841 

15 

7,9905 

14 

7.0000 

0 

4.2019 

10 

1.5980 

10 

2.5536 

20 

8.0000 

0 

2.5148 

13 

6.3240 

26 

9.0000 

0 

1.0947 

17 

1.1985 

34 

1.0000 

1 

1.3130 

27 

1.7240 

42 

Table  6.2(a).  The  mean  and  standard  deviation  of  the  upward  first  passa£e 
time  T„„ 

- OX  .^0 


X 


10 


25 


40 


80 


1.0000 

-1 

1.2307 

-1 

1.2024 

-1 

4.3176 

-1 

4.0435 

-1 

2.0000 

-1 

2.3012 

-1 

2. 3055 

-1 

5.5216 

-1 

5.4151 

-1 

3.0000 

-1 

3.3060 

-1 

3.3343 

-1 

6.3979 

-1 

6.3626 

-1 

4.0000 

-1 

4.2497 

-1 

4.2957 

-1 

7.0793 

-1 

7.0760 

-1 

5.0000 

-1 

5.1379 

-1 

5.1953 

-1 

7.6292 

-1 

7.6339 

-1 

6.0000 

-1 

5.9738 

-1 

6.0393 

-1 

8.0794 

-1 

8.0825 

-1 

7.0000 

-1 

6.7609 

-1 

6.8331 

-1 

8.4513 

-1 

8.4504 

-1 

8.0000 

-1 

7.5036 

-1 

7.5816 

-1 

8.7622 

-1 

8.7567 

-1 

9.0000 

-1 

8.2058 

-1 

8-2891 

-1 

9.0247 

-1 

9.0151 

-1 

1.0000 

0 

8.8710 

-1 

8.9591 

-1 

9.2479 

-1 

9.2351 

-1 

1.2000 

0 

1.0103 

0 

1.0200 

0 

9.6055 

-1 

9.5875 

-1 

1.3000 

0 

1.0675 

0 

1.0776 

0 

9.7500 

-1 

9.7297 

-1 

1.6000 

0 

1,2245 

0 

1.2354 

0 

1.0088 

0 

1.0061 

0 

2.0000 

0 

1.4054 

0 

1.4172 

0 

1.0388 

0 

1.0354 

0 

2.2000 

0 

1.4861 

0 

1.4982 

0 

1.0496 

0 

1.0459 

0 

2.6000 

0 

1.6319 

0 

1.6445 

c 

1.0660 

0 

1.0618 

0 

3.0000 

0 

1.7605 

0 

1.7735 

0 

1.0775 

0 

1.0730 

0 

3.5000 

0 

1.9023 

0 

1.9157 

0 

1.0875 

0 

1.0827 

0 

4.0000 

0 

2.0274 

0 

2.0410 

0 

1.0946 

0 

1.0895 

0 

5.0000 

0 

2.2541 

0 

1.0980 

0 

6.0000 

0 

2.4307 

0 

1.1029 

0 

7.0000 

0 

2.5813 

0 

1.1060 

0 

8.0000 

0 

9.0000 

0 

1.0000 

1 

Table  6.2(b).  The  mean  and  standard 


1.2012  -1 

1.2025  -1 

1.2053  -1 

4.0050  -1 

3.9896  -1 

3.9887  -1 

2.3108  -1 

2.3163  -1 

2.3221  -1 

5.4076  -1 

5.4080  -1 

5.4099  -1 

3.3440  -1 

3.3521  -1 

3.3603  -1 

6.3644  >1 

6.3661  -1 

6.3675  -1 

4.3080  -1 

4.3182  -1 

4.3284  -1 

7.0778  -1 

7.0788  -1 

7.0796  -1 

5.2097  -1 

5.2218  -1 

5.2339  -1 

7.6347  -1 

7.6349  -1 

7.6351  -1 

6.0556  -1 

6.0694  -1 

6.0831  -1 

8.0821  -1 

3.0817  -1 

8.0812  -1 

6.851?  “1 

6.8664  -1 

6.8817  -1 

8.4491  -1 

8.4480  -1 

8.4468  -1 

7.6014  -1 

7.6179  -1 

7.6344  -1 

8.7547  -1 

8.7530  -1 

8.7511  -1 

8.3102  -1 

8.3279  -1 

8.3457  -1 

9.0123  -1 

9.0099  -1 

9.0074  -1 

8.9815  -1 

9.0003  -1 

9.0191  -1 

9.2316  -1 

9.2286  -1 

9.2254  -1 

1.0224  0 

1.0245  0 

1.0266  0 

9.5826  -1 

9.5784  -1 

9.5741  -1 

1.0801  0 

1.0823  0 

1.0844  0 

9.7242  -1 

9.7195  -1 

9.7146  -1  1 

1.2382  0 

1.2405  0 

1.2428  0 

1.0054  0 

1.0047  0 

1.0041  0 

1.4202  0 

1.4227  0 

1.4252  0 

1.0345  0 

1.0337  0 

1.0329  0  I 

1.5013  0 

1.5039  0 

1.5065  0 

1.0450  0 

1.0441  0 

1.0433  0 

1.6477  0 

1.6504  0 

1.6531  0  i 

1.0607  0 

1.0598  0 

1.0589  0 

1.7768  0 

1.7796  0 

1.7824  0 

1.0718  0 

1.0708  0 

1.0697  0  1 

1.9191  0 

1.9219  0 

1.9247  0 

1.0814  0 

1.0803  0 

1.0792  0 

2.0445  0 

2.0474  0 

2.0503  0 

1.0881  0 

1.0870  0 

1.0858  0 

2. 2576  0 

2.2606  0 

2.2635  0 

1.0966  0 

1.0954  0 

1.0941  0 

2.4343  0 

2.4373  0 

2.4403  0 

1.1015  0 

1.1002  0 

1.0989  0 

2.5849  0 

2.5880  0 

2.5910  0 

1.1045  0 

1.1032  0 

1.1019  0 

2.7162  0 

2.7192  0 

2.7223  0  i 

1.1066  0 

1.1052  0 

1.1039  0  I 

2.8354  0 

2.8385  0 

1.1067  0 

1.1053  0  1 

2.9396  0 

2.9427  0  i 

1.1077  0 

1.1063  0 

\ 

eviation  of 

the  downward 

first  passage 

time  T 


■XO 


:ii 


X 

10 

25 

40 

i 

0.10 

0.0000 

0.0000 

0.0000 

0.0000 

0.20 

0.0000 

0.0000 

0.0000 

0.0000 

[ 

0.30 

0.0000 

0.0000 

0.0000 

0.0000 

' 

0.4S 

0.0002 

0.0002 

0.0002 

0.0002 

' 

0.60 

0.0007 

0.0006 

0.0006 

0.0006 

\ 

¥ 

0.75 

0.0014 

0.0014 

0.0014 

0.0014 

t 

0.90 

0.0023 

0.0024 

0.0025 

0.0025 

1  ■ 

1.00 

0.0030 

0.0032 

0.0033 

0.0033 

1 

1.20 

0.0046 

0.0050 

0.0051 

0.0052 

1.40 

0.0064 

0.0070 

0.0071 

0.0073 

1 

t 

1.60 

0.0083 

0.0091 

0.0093 

0.0095 

2.00 

0.0121 

0.0133 

0.0136 

0.0141  j 

2.20 

0.0141 

0.0155 

0.0158 

0.0163  I 

k 

2.60 

0.0180 

0.0198 

0.0203 

0.0206  1 

0.0255  ! 

i 

3.00 

0.0219 

0.0241 

0.0247 

3.40 

0.0258 

0.0284 

0.0291 

0.0300  ! 

3.80 

0.0297 

0.0327 

0.0335 

0.0345  1 

t  ■ 

5.00 

0.0413 

0.0454 

0.0464 

0.0479 

*■ 

Table  6 

3(a).  Sq^(x) 

jL-P-1Io3^ 

XyCt) 

X 

10 

25 

40 

X 

o 

c: 

0.1000 

0.0001 

0.0000 

0.0000 

0.0000 

0.2000 

0.0004 

0.0003 

0.0003 

0.0002 

C.3000 

0.0014 

0.0013 

0.0013 

0.0012 

0.4500 

0.0037 

0.0040 

0.0040 

0.0041 

0.6000 

0.0066 

0.0073 

0.0075 

0.0077 

0.7500 

0.0097 

0.0108 

0.0111 

0.0115 

0.9000 

0.0126 

0.0141 

0.0145 

0.0151 

1.0000 

0.0145 

0.0162 

0.0167 

0.0174 

1.2000 

0.0180 

0.0202 

0.0207 

0.0216 

1.4000 

0.0212 

0.0238 

0.0244 

0.0254 

1.6000 

0.0242 

0.0271 

0.0278 

0.0289 

2.0000 

0.0296 

0.0330 

0.0339 

0.0351 

2.2000 

0.0320 

0.0357 

0.0366 

0.0380 

2.6000 

0.0367 

0.0408 

0.0418 

0.0433 

3.0000 

0.0410 

0.0455 

0.0467 

0.0483 

3.4000 

0.0451 

0.0500 

0.0513 

0.0531 

3.8000 

0.0491 

0.0544 

0.0558 

0.0577 

5.0000 

0.0607 

0.0671 

0.0687 

0.0710 

Table  6.3(b).  Sjj(x)  =  PlTj^  £  x] 
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L 


X 

10 

25 

40 

0.10 

0.0080 

0.0C81 

0.0076 

0.0067 

0.20 

0.0242 

0.0272 

0.0279 

0.0288 

0.30 

0.0389 

0.0445 

0.0459 

0.0479 

0.45 

0.0557 

0.0636 

0.0656 

0.0684 

0.60 

0.0675 

0.0767 

0.0790 

0.0822 

0.75 

0.0763 

0.0863 

0.0888 

0.0922 

0.90 

0.0830 

0.0935 

0.0962 

0.0996 

1.00 

0.0866 

0.0975 

0.1002 

0.1040 

1.20 

0.0927 

0.1040 

0.1069 

0.1108 

1.40 

0.0975 

0.1092 

0.1122 

0.1162 

1.60 

0.1015 

0.1135 

0.1166 

0.1207 

2.00 

0.1080 

0.1206 

0.1237 

0.1280 

2.20 

0.1108 

0.1235 

0.1268 

0.1311 

2.60 

0.1157 

0.1288 

0.1322 

0.1367 

3.00 

0.1200 

0.1336 

0.1370 

0.1416 

3.40 

0.1241 

0.1379 

0.1414 

0.1462 

3.80 

0.1279 

0.1421 

0.1456 

0.1505 

5.00 

0.1387 

0.1538 

0.1576 

0.1628 

Table  6.3(c). 


i23iiiL 


X 

10 

25 

40 

o 

X 

0.10 

0.0001 

0.0000 

0.0000 

0.0000 

0.20 

0.0019 

0.0011 

0.0009 

0.0006 

0.30 

0.0062 

0.0048 

0.0044 

0.0038 

0.45 

0.0170 

0.0152 

0.0147 

0.0140 

0.60 

0.0309 

0.0290 

0.0286 

0.0278 

0.75 

0.0461 

0.0443 

0.0438 

0.0431 

0.90 

0.0617 

0.0599 

0.0594 

0.0588 

1.00 

0.0720 

0.0702 

0.0698 

0.0691 

1.20 

0.0923 

0.0905 

0.0901 

0.0894 

1.40 

0.1119 

0.1100 

0.1096 

0.1089 

1.60 

0.1308 

0.1288 

0.1284 

0.1276 

2.00 

0.1665 

0.1644 

0.1639 

0.1631 

2.20 

0.1835 

0.1813 

0.1808 

0.1800 

2.60 

0.2161 

0.2137 

0.2132 

0.2123 

3.00 

0.2470 

0.2445 

0.2439 

0.2430 

3.40 

0.2764 

0.2738 

0.2732 

0.2722 

3.80 

0.3046 

0.3018 

0.3012 

0.3002 

5.00 

0.3824 

0.3793 

0.3786 

0.3775 

Table  6.3Cd). 


§02^^ 


X 

10 

25 

40 

0.10 

0.0203 

0.0160 

0.0146 

0.0119 

0.20 

0.0591 

0.0559 

0.0551 

0.0538 

0.30 

0.0955 

0.0940 

0.0938 

0.0934 

0.4S 

0.1399 

0.1397 

0.1397 

0.1398 

0.60 

0.1744 

0.1746 

0.1747 

0.1749 

0.75 

0.2021 

0.2024 

0.2025 

0.2027 

0.90 

0.2253 

0.2256 

0.2256 

0.2258 

1.00 

0.2389 

0.2391 

0.2392 

0.2393 

1.20 

0.2628 

0.2630 

0.2630 

0.2631 

1.40 

0.2837 

0.2837 

0.2837 

0.2837 

1.60 

0.3023 

0.3022 

0.3022 

0.3022 

2.00 

0.3351 

0.3348 

0.3347 

0.3346 

2.20 

0.3499 

0.3495 

0.3494 

0.3493 

2.60 

0.3775 

0.3768 

0.3767 

0.3765 

3.00 

0.4029 

0.4021 

0.4019 

0.4017 

3.40 

0.4268 

0.4258 

0.4256 

0.4253 

3.80 

0.4494 

0.4483 

0.4480 

0.4477 

5.00 

0.5113 

0..5098 

0.5095 

0.5090 

Table  6 

.3(e).  Sj.,(x)  =  P[Tj., 

X 

10 

25 

40 

w  1 

o 

X 

0.10 

0.0386 

0.0285 

0.0257 

0.0204 

0.20 

0.1103 

0.1005 

0.0981 

0.0941 

0.30 

0.1787 

0.1710 

0.1692 

0.1663 

0.45 

0.2640 

0.2579 

0.2564 

0.2539 

0.60 

0.3319 

0.3263 

0.3249 

0.3226 

0.75 

0.3876 

0.3821 

0.3807 

0.3785 

0.90 

0.4346 

0.4292 

0.4278 

0.4256 

1.00 

0.4623 

0.4569 

0.4556 

0.4534 

1.20 

0.5113 

0.5060 

0.5047 

0.5025 

1.40 

0.5536 

0.5483 

0.5470 

0.5449 

1.60 

0.5908 

0.5857 

0.5844 

0.5823 

2.00 

0.6540 

0.6491 

0.6479 

0.6459 

2.20 

0.6812 

0.6765 

0.6753 

0.6733 

2.60 

0.7287 

0.7243 

0.7232 

0.7213 

3.00 

0.7687 

0.7646 

0.7636 

0.7619 

3.40 

0.6027 

0.7989 

0.7979 

0.7963 

3.80 

0.8315 

0.8281 

0.8272 

0.8257 

5.00 

0.8951 

0.8925 

0.8918 

0.8907 

Table  6.3(fJ.  Sq^(x)  ^  PfTp^  s  x] 


34 


0 


XyCt) 


X 

10 

25 

40 

0.1000 

0.1294 

0.1028 

0.0951 

0.0808 

0.2000 

0.3501 

0.3406 

0.3389 

0.3367 

0.3000 

0.5384 

0.5413 

0.5424 

0.5447 

O.HSQO 

0.7308 

0.7368 

0.7383 

0.7408 

0.6000 

0.8436 

0.8477 

0.8487 

0.8503 

0.7500 

0.9068 

0.9111 

0.9116 

0.9125 

0.9000 

0.9466 

0.9478 

0.9480 

0.9484 

1.0000 

0.9626 

0.9633 

0.9634 

0.9637 

1.2000 

0.9816 

0.9818 

0.9818 

0.9819 

1.4000 

0.9909 

0.9909 

0.9909 

0.9909 

1.6000 

0.9955 

0.9955 

0.9954 

0.9940 

2.0000 

0.9989 

0.9989 

0.9989 

0.9988 

2.2000 

0.9995 

0.9994 

0.9994 

0.9999 

2.6000 

0.9999 

0.9999 

0.9999 

0.9999 

3.0000 

1.0000 

1.0000 

1.0000 

1.0000 

3.4000 

1.0000 

1.0000 

1.0000 

1.0000 

3.8000 

1.0000 

1.0000 

1.0000 

1.0000 

5.0000 

1.0000 

1.0000 

1.0000 

1.0000 

Table  6. 

2(x)  =  P[T_ 

3.-2  “ 

Xv(t) 


X 

1C 

25 

40 

o 

X 

0.10 

0.0025 

0.0007 

0.0004 

0.0001 

0.20 

0'.0268 

0.0173 

C.0148 

0.0103 

0.30 

0.0857 

0.0718 

0.0680 

0.0612 

0.45 

0.2196 

0.2086 

0.2058 

0.2011 

0.60 

0.3689 

0.3635 

0.3622 

0.3602 

0.75 

0.5060 

0.5045 

0.5042 

0.5038 

0.90 

0.6209 

0.6214 

0.6216 

0.6219 

1.00 

0.6845 

0.6856 

0.6859 

0.6864 

1.20 

0.7839 

0.7852 

0.7856 

0.7862 

1.40 

0.8532 

0.8545 

0.8548 

0.8553 

1.60 

0.9009 

0.9018 

0.9020 

0.9024 

2.00 

0.9551 

0.9556 

0.9557 

0.9559 

2.20 

0.9699 

0.9702 

0.9703 

0.9704 

2.60 

0.9864 

0.9866 

0.9866 

0.9867 

3.00 

0.9939 

0.9940 

0.9940 

0.9940 

3.40 

0.9973 

0.9973 

0.9973 

0.9973 

3.80 

0.9986 

0.9988 

0.9988 

0.9988 

5.00 

0.9999 

0.9999 

0.9999 

0.9999 

Table  6. 

,3(h).  S_3_ 

<  X] 
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t-. : 


I  ' 


r 


\Ct) 


X 

10 

25 

40 

X 

o 

c: 

rt 

0.10 

0.0012 

0.0003 

0.0002 

0.0000 

0.20 

0.0133 

0.0079 

0. 0056 

0.0043 

0.30 

0.0433 

0.0341 

''.0316 

0.0274 

0.45 

0.1160 

0.1053 

1025 

0.0978 

0.60 

0.2051 

0.1956 

0.1932 

0.1892 

0.75 

0.2966 

0.2890 

0.2871 

0.2838 

0.90 

0.3840 

0.377? 

0.3761 

0.3735 

1.00 

0.4379 

0.4325 

0.4311 

0.4288 

1.20 

0.5346 

0.5304 

0.5293 

0.5276 

1.40 

0.6163 

0.6130 

0.6122 

0.6106 

1.60 

0.6844 

0.6819 

0.6612 

0.6801 

2.00 

0.7875 

0.7858 

-  0.7854 

0.7847 

2.20 

0.8258 

0.8245 

0.8241 

0.8236 

2.60 

0.8831 

0.8822 

0.8820 

0.8816 

3.00 

0-9216 

0.9210. 

0.9208 

0.9206 

3.40 

0.9474 

0.9470 

0.9469 

0.9468 

3.80 

0.9647 

0.9645 

0.9644 

0.9643 

5.00 

0.9894 

0.9693 

0.9893 

0.9892 

Table 

6.3(i).  S_2 

<  X] 

1 

1 

1 

X 

< 

rt 

! 

i 

i 

i 

1 

X 

10 

25 

40 

X 

o 

c: 

ft 

0.10 

0.0629 

0.0465 

0.0420 

0.0336 

0.20 

0.1773 

0.1632 

0.1597 

0.1539 

0.30 

0.2853 

0.2759 

0.2737 

0.2701 

0.45 

0.4171 

0.4115 

0.4101 

0.4078 

0.60 

0.5183 

0.5142 

0.5132 

0.5115 

0.75 

0.5973 

0.5941 

0.5933 

0.5920 

0.90 

0.6607 

0.6581 

0.6574 

0.6563 

1.00 

0.6963 

0.6939 

0.6934 

0.6924 

1.20 

0.7552 

0.7534 

0.7529 

0.7521 

1.40 

0.8016 

0.8001 

0.7998 

0.7991 

1.60 

0.8387 

0.8375 

0.8372 

0.8367 

2.00 

0.8927 

0.8919 

0.8917 

0.8914 

2.20 

0.9123 

0.9117 

0.9115 

0.9112 

2.60 

0.9414 

0.9409 

0.9408 

0.9406 

3.00 

0.9607 

0.9604 

0.9604 

0.9602 

3.40 

0.9737 

0.9735 

0.9734 

0.9734 

3.80 

0.9824 

0.9822 

0.9822 

0.9821 

5.00 

0.9947 

0.9947 

0.9946 

0.9946 

Table 

Table  6.3(j).  S_j^q(x)  "  ^[^-1,0^  x] 


i 


I 

A 


36 


Xy(t) 


X 

10 

25 

4(! 

0.10 

0.0000 

0.0000 

0.0000 

0.0000 

0.20 

0.0005 

0.0002 

0.0001 

0.0000 

0.30 

0.0042 

0.0021 

0.0017 

0.0009 

0.45 

0.0233 

0.0172 

0.0157 

0.0130 

0.60 

0.0640 

0.0552 

0.0529 

0.0489 

0.75 

0.1243 

0.1146 

0.1121 

0.1079 

0.90 

0.1972 

0.1881 

0.1858 

0.1819 

1.00 

0.2494 

0.2410 

0.2389 

0.2353 

1.20 

0.3547 

0.3479 

0.3462 

0.3433 

1.40 

0.4541 

0.4488 

0.4475 

0.4452 

1.60 

0.5432 

0.5391 

0.5380 

0.5363 

2.00 

0.6863 

0.6837 

0.6831 

0.6820 

2.20 

0.7415 

0.7394 

0.7389 

0.7380 

2.60 

0.8255 

0.8241  • 

0.8238 

0.8232 

3.00 

0.9826 

0.8817 

0.8815 

0.8611 

3.40 

0.9212 

0.9206 

0.9205 

0.9202 

3.80 

0.9471 

0.9467 

0.9466 

0.9465 

5.00 

0.9841 

0.9840 

0.9839 

0.9839 

Table  6.3(k).  S_2  ^(x)  =  P[T_3  ^  s  x] 


X 

10 

25 

40 

X 

O 

c: 

r* 

0.10 

0.0932 

0.0711 

0.0648 

0.0531 

0.20  • 

0.2582 

0.2440 

0.2406 

0.2353 

0.30 

0.4084 

0.4032 

0.4021 

0.4006 

0.45 

0.5801 

0.5802 

0.5803 

0.5806 

0.60 

0.6992 

0.7006 

0.7010 

0.7016 

0.75 

0.7825 

0.7840 

0.7844 

0.7850 

0.90 

0.8416 

0.8429 

0.8432 

0.8438 

1.00 

0.8713 

0.8725 

0.8728 

0.8733 

1.20 

0.9147 

0.9156 

0.9158 

0.9161 

1.40 

0.9432 

0.9438 

0.9440 

0.9442 

1.60 

0.9621 

0.9625 

0.9626 

0.9628 

2.00 

0.9831 

0.9832 

0.9833 

0.9834 

2.20 

0.9887 

0.9888 

0.9888 

0.9889 

2.60 

0.9949 

0.9950 

0.9950 

0.9950 

3.00 

0.9977 

0.9977 

0.9977 

0.9978 

3.40 

0.9930 

0.9990 

0.9990 

0.9990 

3.80 

0.9995 

0.9995 

0.9995 

0.9995 

5.00 

1.0000 

1.0000 

1.0000 

1.0000 

Table6.3(£).  S_2  _ 

j(x)  =  PtT.2, 

-1  ^ 

